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‘@} HISTORY ; Detailed Current Edit History 
(5) PECL ARAT 1ONS 

(4) MTHSCGSIN = G COMPLEX*1 SING 

(5) 1 $ MTHSCGCOS = G COMPLEX*16 Cosine 

(6) 1 WORKER = do all the work 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


a? fe ucche a gle UNDER A LICENSE AND MAY BE USED AND COPIED 


® 
® 
® 
® 
® 
® 
ACCORDANC OF SUCH LICE » 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
TRANSFERRED. : 
® 
THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
® 
* 
® 
& 
® 
* 
® 
® 


AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
CORPORATION. 


DIGITAL ASSUMES NO RESPONS sia 


IBILIT 
SOFTWARE ON EQUIPMENT WHICH 


FOR THE USE OR emma OF ITS 


LITY 
S NOT SUPPLIED BY DIGITAL 
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fACILITY: MATH LIBRARY 


“ABSTRACT: 
Return the SINE of a G COMPLEX*16 number 
Return the COSINE of a G COMPLEX*16 number 


VERSION: 1 
HISTORY: 


AUTHOR: 
Steven B. Lionel, 26-July-1979 


MODIFIED BY: 


2 ants cman 
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M005 6 SPRY ; Detailed Current Edit History g-3 sets 7 9:99:98 ERTHRTL « SRCIJMTHCGSINC.MAR; 1 ~__ 3) 
-SBTTL HISTORY 3; Detailed Current Edit History 
oi 3 Edit History 
Bh gi Y 1-00) = Adapted from MTHSCSINCOS version 1-002. SBL 26-July-1979 
000 3 O02 = Use MTHSGEXP_R6. SBL 14-Dec-1979 


g Ss 
MTHSCGSINCOS G COMPLEX*16 Sine and Cosine 16-SEP-1984 01:09:4 AX/VMS Macro V04-00 Page 
pa + SEE ARAT TONG paises o: 7369 YRTHRTL. SREIRTHCGSINC MAR: 1 ° 
23 »SBTTL DECLARATIONS 
61; 
86 3; INCLUDE FILES: 
65 ; 
64 
3 65 ; 
0 68 3; EXTERNAL SYMBOLS: 
00 6 -DSABL G 
00 68 eEXTRN MTHSGSIN_R7 
4 8 eEXTRN MTHSGCOS_R7 
B88 0 sEXTRN MTHSGEXP_R6 
00 71 
0000 ue 3 
0000 73 ; MACROS: 
444 74 ; 
000 75 
0000 6 ; 
0000 77 ; PSECT DECLARATIONS: 
vest °° ee -PSECT _MTHSCODE PIC, SHR, LONG, EXE, NOWRT 
$000 80 ; 
0000 81 ; EQUATED SYMBOLS: 
0000 Ht 
0000 83 ; 
0000 84 ; OWN STORAGE: 
0000 85 ;: NONE 
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G COMPLEX 
MTHSCGSIN = G COMPLEX®*16 SINE 6-SEP-1984 MTMRTL. SRCIMTHCGSINC MAR: 1 (4) Tab 
i -SBTTL MTHSCGSIN = G COMPLEX*16 SINE | 
r++ 
a 3 FUNCTIONAL DESCRIPTION: 
0 3 : MTHSCGSIN computes the SINE of a G COMPLEX*16 number (r, i) as | 
Be s result = (SIN(r) * COSH(i), COS(r) * SINH(i)) | 
4 38 3; CALLING SEQUENCE: 
38 $3 : CALL MTHSCGSIN (result.wgc.r, arg.rgc.r) | 
000 100 ; INPUT PARAMETERS: 
00000008 448 19} arg = 8 ; G COMPLEX*16 argument by reference 
0000 198 3; IMPLICIT INPUTS: 
44 : NONE 
000 105; 
4 138 ; OUTPUT PARAMETERS: 
00000004 448 198 result = 4 ; G COMPLEX*16 result by reference 
0000 110 ; IMPLICIT OUTPUTS: 
0000 111; NONE 
0000 WN 3 
0000 113 ; COMPLETION CODES: 
0000 114; NONE 
0000 115; 
0000 116 ; SIDE EFFECTS: = 
0000 117; Signals: Reserved Operand if r or i are invalid (-0.0) 
0000 118; or SINSIGLOS if iri > 2*PI*2**31, 
0000 119; Floating Overflow if i > 88.028. 
0000 19 ; 
0000 121 ;-- 
0000 1$ 
0000 1$ 
OOF C 944 124 ENTRY MTHSCGSIN, “M<R2,R3,R4,R5,R6,R7> 
OOOO0O3E'EF 16 008 125 JS WORKER : RO-R1 =SIN(r) 
000 126 : R2-R3 = COS(r) 
0008 127 : R4-R5 = SINH(i) 
0008 128 : R6-R7 = COSH(i) 
50 56 44FD 000 129 ; RO-R1 = SIN(r) * COSH(i) 
52. 54 44FD st 130 3: R2-R3 = COS(r) * SINH(i) 
54 04 AC 34 1 131 3: get result address 
84 50 D 0014 1 $ 3 
64 52 7D Q017 1 ; Store taee (nar? part 
04 OO1A 134 


| 
tore real part 


—— ——s a 
| 
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M4 - LEX®* osine -SEP- ré1: MTHRTL.SRCIMTHCGSINC.MAR; (5) 
M005 aThscethe ~ G6 COMPLEX#{o Cost Breer bbe S1ist;09 EMTARTE eaeaarie gr 
18 1 § -SBTTL MTHSCGCOS - G COMPLEX*16 Cosine 
1B 138 :4s 
18 133 : FUNCTIONAL DESCRIPTION: 
18 13) ; MTHSCGCOS computes the COSINE of G COMPLEX*16 number (r, i) as follows: 
i8 148 ; result = (COS(r) * COSH(1), -SIN(r) * SINK(-i)) 
18 145 : CALLING SEQUENCE: 
i8 178 3 CALL MTHSCGCOS (result.wgc.r, arg.rgc.r) 
1B 148; 
018 128 ; INPUT PARAMETERS: 
00000008 aie 139 arg = 8 ; G COMPLEX*16 argument by reference 
018 136 : IMPLICIT INPUTS: 
Soe 18a as 
0018 135 > OUTPUT PARAMETERS: 
00000004 Sole 139 , result = 4 ; G COMPLEX*16 result by reference 
0018 159 : IMPLICIT OUTPUTS: 
os ie; 
0018 168 > COMPLETION CODES: 
OO1k 164 te: 
0018 165 : SIDE EFFECTS: 
018 166; Signals: Reserved Operand if r or i are invalid (-0.0) 
0018 167 ; MTH$_SINSIGLOS if iri > 2*PI #22231, 
001B 168; Floating Overflow if i > 88.028. 
0018 199 3 
018 170 ;-- 
Bie 
OOF C 018 153 -ENTRY MTHSCGCOS, “M<R2,R3,R4,R5,R6,R7> 
QOOOOO3E'EF 16 91 174 JSB WORKER : RO-Ri = SIN(r) 
Hy = = r 
80 176 3 R4-R5 = SINH(i) 
177 3 R6-R7 = COSH(i) 
D 00 178 MULG2 R2, R6 3 R6-R7 = COS(r) * COSH(i) 
D 17 MNEG Re. RO 3; RO-R1 = -SIN(r) 
D 180 MULG RO, R4, R2 3 R2-R3 = -SIN(r) * SINH(i) 
D ; 181 MOVQ R6, R : RO-R1 = COS(r) * COSH(i) 
0 1 ¢ MOVL result(AP), R4 3 get result address 
D . 1 MOva Ro. (R4)+ 3; Store real part 
? . : : peye R2, (R4) ; Store imaginary part 
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M005 WORKER = do all the work -SEP-1984 ot 7:9 Ya ANTe SAE IRTHCCS GHC MAR: 1 ~— (8) 1-0) 
—E 187 ~SBTTL WORKER = do all the work 
ie. 
ce. : Setup error handler 
E 194 3 Eonpute: 
E 138 ; RO-R1 = SIN(r) 
—E 193; R2-R3 = COS(r) 
E 194; R4-R5 = SINH(i) 
—E 195; R6-R7 = COSH(i) 
E 199° 
8 E 138 WORKER: 
8 : 19 MTHSFLAG_JACKET ; set up error handler 
6D 00000000'GF 9E 00 MOVAB G“MTHSSJACKET_HND, (FP) 
1 ; set handler address to jacket 
bee 3; handler 
50 08 AC dO 045 00 MOVL (AP), RO 3; RO => ¢ i) 
$0 08 AO SOFD 0049 501 MOVE (RO) SORT of 
90000000" EF 16 004E 0¢ JSB MTHSGEXP_R6 > RO-R1 = EXP(i) 
5 08 0 47FD 424 0% DIVG3 RO, #1, R2 3 R2-R3 = EXP(-i) 
54 50 52 43FD $089 05 SUBG R2, RO, R4 3 R4-R5 = EXPC(I) = EXP(-i) 
7E 54 00 45FD O68 89 MULG #0.5, R4, -(SP) 3; (SP) = SINH(i) 
54 50 52 41FD 06 08 ADDG R2, RO, R4 3 R4-R5 = EXPCi) + EXP(-i) 
7E 54 00 45FD $04 MULG #0.5, R4, -(SP) ; (SP) = COSH(i) 
08 BC SOFD pee 11 MOVG @arg(AP), RO 3; RO-RI = Fr 
sees tall 16 00 § \¢ B MTHSGCOS_R7 3 RO-R1 = COS(r) 
7E 0 7D $3 17 MOVQ RO, -(SPY ; (SP) = COS(r) 
50 08 BC 7D ae 15 MOva @arg(AP), RO 3 RO-RI =r 
OQOOOOOO0'EF 16 7F 1 JSB MTHSGSIN_R7 3 RO-R1 = SIN(r) 
52 8E 7D 85 1 MOVQ (SP)+, R2 3 R2-R3 = COS(r) 
56 8E 7D 088 18 MOVa (SP)+, R6 3: R6-R7 = COSH(i) 
54 8E 7D bbe i? MOVa (SP)+, R4 3 R4-R5 = SINH(i) 
OO8E 1 
09 F § 
O8F 
+H 4 
008F 5 
| 
| 
J 
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ARG = 00000008 
MIHSSJACKET_HND eeeeeere 


Cross-reference output 
Assembler run totals 


1 
MTHSCGCOS 1B RG 1 
MTHSCGSIN Mints RG 1 
MTHSGCOS_R7 eeeerere x 6 
MTHSGEXP_R eteeeree = § X 
MTHSGSIN_R eeeerere xX 
RESULT = 44 
WORKER OOOOOSE R 01 
t+aeeeceaer ene ceaeaca} 
! Psect synopsis ! 
pecan www nm weer enen} 
PSECT name Allocation PSECT No. Attributes 
- ABS. 444444 ( 9.3 00 « 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE QOOCOOSF (« 1435.) O01¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
grew r er ere ne seco ewe em ee om Se 
; Performance indicators ; 
Phase Page faults CPU Time Elapsed Time 
Initialization 37 0: 8: 0.11 8:00:00.6 
Command processing 119 0:00:00 4 0:00:04.4 
Pass 1 91 0:00:00.7 :00:04.4 
Symbol table sort 0 :00: 3 00 :00:00.00 
Symbol table output 3 :00 G3 ; 6; 0:09 
able outpu :00:00. 
f 6) 00;00:00.04 
3 00: 0:17.72 


ooooceo 
ooooocoeo 


0 
00: 
Ba: : 
Psect synopsis output 3: : Q 
30 00:00; 


The working set Limit was 750 pages. 

3610 bytes (8 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 9 non-local and 0 local symbols. 
source Lines were read in Pass 1, producing 14 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 


peseecee ee ne eco er acne sees 


Macro Library name Macros defined 
“S$255SDUAZB:(SYSLIBJSTARLET.MLBs20 tS 0 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LISS:MTHCGSINC/OBJ=OBJ$:MTHCGSINC MSRCS$:MTHJACKE T/UPDATE=(ENHS:MTHJACKET) +MS 
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